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Results are presented from a Monte Carlo simulation of the Nambu - Jona-Lasinio model with SU(2)(8)SU(2) 
chiral symmetry and Nf = 2 flavors of fermion. We show that fits to the equation of state are sensitive to the 
shape and extent of the assumed scaling region; the best flts favour a triviality scenario predicted by the large- A^/ 
approximation, in which the scalar degrees of freedom are fermion - anti-fermion composites with wavefunction 
renormalisation constant vanishing logarithmically in the continuum limit. 
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In this talk I report on work done with John 
Kogut; a more detailed presentation has ap- 
peared . We have studied the critical behaviour 
of the four-dimensional Nambu - Jona-Lasinio 
(NJL) model, whose continuum Lagrangian den- 
sity is given by 
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where i runs over Nj distinct fermion flavors, 
and in the second line we introduce an auxil- 
iary scalar field $ = cr + in.T proportional to 
an SU(2) matrix. The Lagrangian (]l|) has an 
SU(2)^ (g) SU(2)^ chiral symmetry i/'l i-^ UtpL, 
tfjji 1-^ VipR, $ 1-^ V^W, which is broken explic- 
ity by the bare fermion mass, and spontaneously 
by the development of a condensate E = (cr) ^ 0. 

The NJL model has been proposed to describe 
both strong |^ and weak js) interaction physics; 
in the latter case the Higgs is a composite // 
state. Both the NJL model and the conventional 
model of the Higgs sector, namely the 0(4) linear 
sigma model, are believed trivial in the sense that 
neither has an interacting continuum limit. It has 
been shown that the lattice-regularised models 
have equivalent predictive power for electroweak 
phenomenology Q. 

In this study we aim to characterise the trivi- 
ality of the NJL model by fitting numerical data 



for S to the model equation of state (EOS) 
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The logarithms are corrections to the EOS pre- 
dicted by mean field theory. Their presence sig- 
nals triviality: the interacting theory becomes ill- 
defined in the continuum limit S ^ (in lat- 
tice units). Eq. (|^) has wide applicability; for 
the O(A^) sigma model the exponents qi^2 can be 
calculated using RG-improved perturbation the- 
ory l^j. For the NJL model, they can be estimated 
to leading order in the expansion ||], and 

are qualitatively distinct (see table). We wish 
to understand whether this difference is real or 
merely an artifact of the approximation. 

We have simulated a discretised version of (|^) 



on a 16 lattice for N 



f 



2 111. Since we use 



staggered fermions, this necessitates a fractional 
number of lattice flavors, and hence the use of a 
hybrid molecular dynamics algorithm. We have 
accumulated a high statistics dataset with bare 
mass values in the range m = 0.05, . . . , 0.0025. 
We also performed simulations directly in the chi- 
ral limit m = 0. In this case S vanishes in a finite 
volume, so instead we measured the mean field 
|$| defined via 
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In order to find a satisfactory fit to (g) it is 
necessary to make some assumptions about the 
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Table 1 

Values for the correction exponents 





Analytical predictions 




Numerical fits 






0(4) sigma model 


large- iV/ NJL 


ferromagnetic 


composite 




0.5 


0.0 


0.786(19) 


0.016(11) 




1.0 


-1.0 


0.372(17) 


-0.553(18) 
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Figure 1. Fisher plot of the ferromagnetic fit. Figure 2. Fisher plot of the composite fit. 



size and shape of the scaling window about the 
critical point, and hence which data to include. 
We have found that our results for qi depend very 
sensitively on the window chosen. For instance, if 
we include all mass values, but exclude couplings 
outside the range [0.52,0.55] then we obtain val- 
ues, shown in the table, which qualitatively re- 
semble the analytic prediction for the 0(4) sigma 
model: we label this fit "ferromagnetic". If, on 
the other hand, we include all coupling values 
but exclude masses greater than 0.01, then the 
fit changes, the value of qi becoming almost zero, 
and q2 becoming negative. We refer to this fit 
as "composite" because it agrees with the 1/iV/ 
treatment of the NJL model. The results of the 
fits are displayed in Figs. |^, || in the form of a 



Fisher plot of vs. m/S. This would yield 
straight lines of constant coupling as to — > for 
a mean field EOS. 

It would appear that logarithmic corrections 
are hard to pin down from EOS fits alone. Both 
fits have sufficiently large to be questionable. 
Fortunately there is an an alternative probe, in- 
volving data for |<I>| taken at to = 0. On a finite 
system we expect |<I>| to exceed Eq, the true or- 
der parameter obtained by extrapolating the EOS 
to the chiral limit, because \^\ also receives con- 
tributions from fiuctuations of Goldstone modes, 
which average to zero only in the thermodynamic 
limit. The difference A = |$| — Eq can be cal- 
culated using chiral perturbation theory to be of 
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Figure 3. A vs. l/g^: 
fit for 



the form 



the vertical line shows the 
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where Z is the wavefunction renormalisation con- 
stant, defined as the coefficient of l/p^ in the 
pion propagator in the IR limit. The values of 
A obtained using the chiral extrapolation of the 
composite EOS fit obtained by setting qi = 0, 
(72 = — 1 are shown as a function of in 
Fig. H (A extracted using the ferromagnetic EOS 
changes sign over the range and makes no physi- 
cal sense). 

At first sight the data look incompatible with 
(H) , since A is approximately constant for 1 /g^ e 
[0.45,0.52], while Eg falls by a factor of 3 in the 
same range. However, we must take into account 
the dependence of Z on Sq. For models such as 
the sigma model in which the Higgs is elemen- 
tary, Z is perturbatively close to 1 approaching 
the continuum limit, and hence A oc I/Eq. For 
a composite Higgs, however, we expect Z to van- 
ish logarithmically in the same limit Dg, and 
hence A cx 1/Eol InSol- The latter form yields a 
reasonable fit as shown in Fig. |[ 



To conclude, this independent study of the fi- 
nite volume correction seems to prefer the com- 
posite EOS predicted in the large-TV/ approach. 
It remains an interesting open question whether 
this result is generic to four-dimensional models 
in which scalars are composite, or whether other 
universality classes exist defined by alternative 
lattice formulations ||]. 
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